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Abstract: Drought is a common abiotic stress that considerably limits crop production. The objective of 
this study is to explore the influence of water deficiency on the yield, physiologic and metabolomic 
attributes in upland cotton cultivars (Gossypium hirsutum L). Cotton cultivars, 'Ishonch' and "Tashkent-6' 
were selected to study the relationships among their physiologic, metabolomic and yield attributes during 
water deficiency. Deficit irrigation was designed by modifying the traditional watering protocol to reduce 
water use. Results indicate that cotton cultivars respond differently to water deficit stress. Water deficit 
significantly influenced plant height, the number of internodes, and sympodial branches in both cultivars. 
However, yield components such as the number of bolls, boll seed, lint mass, and individual plant yield 
were significantly reduced only in "Tashkent-6'. The leaf area decreased and the specific leaf weight 
increased in 'Ishonch' under deficit irrigation conditions. However, "Tashkent-6' demonstrated significant 
water loss compared to 'Ishonch', and both cultivars showed reduced transpiration rates. Untargeted 


metabolite profiles of leaves showed clear separation in 'Ishonch', but not in "Tashkent-6' under deficit 
irrigation compared to full irrigation. The individual metabolites such as proline and galactinol showed 
strong association with yield under water deficit stress. Moreover, this study indicates that leaf area and 
transpiration intensity influence yield during water deficiency. In summary, the correlation among 
morpho-physiologic, metabolic, and yield components significantly varied between the two cultivars under 
water deficiency. The flowering stage was sensitive to water stress for both cultivars. The direct 
relationship between physiology, metabolism, and yield may be a useful selection criterion for deter mining 
candidate parents for cotton drought tolerance breeding. 
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1 Introduction 


As sessile organisms, plants are continuously exposed to diverse environmental stresses. Of these 
stresses, drought is the most significant because it limits flora diversification and agricultural 
productivity. Drought currently reduces crop yield more than 50% worldwide (Boyer, 1982; 
Bohnert et al., 1995; Baboev et al., 2017). Drought-areas are expanding (Dai, 2011) and 
envisioned to be the worst in the Central and Southwest Asia by 2050 (Agrawala et al., 2001). 
Cotton is considered as 'white gold' and significantly contributes to the economy in cotton 
producing countries. Cotton cultivars differed among each other by their economic valuable traits 
that are main features of both fiber productivity and cooking oil. Additionally, cotton seeds have 
high protein content and are an excellent feed source for livestock. 

Water deficiency stress significantly affects the fiber productivity, quality, and yield traits. 
Uzbekistan is reported to have 4.31x10° hm? of irrigated lands which is about 10% of the total 
lands of the country (Kulmatov, 2014). Water resources management needs to integrate more 
feasible and optimal irrigation methods to control deficit irrigation (AbdelGadir et al., 2012; 
Sampathkumar et al., 2013; Snowden et al., 2013). The deficit irrigation approach is used to 
increase crop water use efficiency by reducing water or the number of irrigation frequencies 
(Kirda, 2002). Considering that, there is a need to drought tolerant cotton varieties to produce 
more seeds and fiber yield per plant without reduction of fiber quality at reduced irrigation 
frequencies as well drought stress conditions. 

Insufficient soil moisture can affect developing organs especially during blooming, flowering, 
and fruiting stages that negatively affect plant morphological traits and yield components 
(McMichael and Hesketh, 1982; Pettigrew, 2004; Soomro et al., 2011). Water deficiency can 
reduce cell and leaf expansion, stem elongation, and leaf area index (Gerik et al., 1996). Plant leaf 
area and reduction of initiation of new leaves are not only attenuated (Pettigrew, 2004), but leaf 
water potential, stomata closure and conductance are also reduced (Bohnert et al., 1995; Bray, 
1997). Seed germination and seedling establishment are the most critical stages in a water-limited 
environment (Saini and Lalonde, 1997) as well as cotton boll development during flowering 
period (Orgaz et al., 1992). 

Drought stress induces the expression of stress-related transcription factors and genes, such 
as reactive oxygen species (ROS), scavenging, abscisic acid (ABA) or mitogen-activated protein 
kinase signaling genes that activate various drought-related pathways to induce tolerance in the 
plant (Ullah et al., 2017). Various genetic, physiologic and biochemical markers have been used 
in breeding programs for drought tolerance of cotton (Pace et al., 1999; Basal et al., 2009). There 
are several studies on genome, epigenome and transcriptome level to identify drought tolerance 
associated quantitative trait loci (QTLs), genes and transcripts in cotton (Chen et al., 2013; Lu et 
al., 2017; Ullah et al., 2017). 

Morphologic and physiologic changes during insufficient water supply lead to genome 
reorganizations and large-scale alterations in transcriptome by suppressing or enhancing genes 
that result in metabolome alterations (Obata and Fernie, 2012; Witt et al., 2012; Pucholt et 
al., 2015). Metabolite markers have been extensively used in screening plants for tolerance and 
responses to various environmental stresses (Witt et al., 2012; Meyer et al., 2014). Metabolites as 
a final product in gene-to-metabolite concept represent alteration in phenotype than genes that 
might be used as predictive markers (Saxena and Cramer, 2013). Witt et al. (2012) demonstrated 
several metabolites that displayed specific responses during drought stress, levels of which were 
correlated well with certain physiological traits. 

Physiologic indices like leaf photosynthesis, stomatal conductance, and transpiration rate 
showed their potential as indicators for drought tolerance in cotton (Laffray and Louguet, 1990; 
Leidi et al., 1993; Nepomuceno et al., 1998; Ullah et al., 2008). Analyses of Gossypium 
barbadense (Siv'on) and Gossypium hirsutum (F-177) also showed higher mineral and metabolite 
content along with greater water use efficiency in F-177 (Levi et al., 2011). Changes in polyamine 
concentrations are common plant response to variety abiotic stresses. Recent studies in 
polyamines (putrescine, spermidine, and spermine) indicated different concentrations in 
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drought-tolerant and drought-sensitive G. hirsutum varieties (Loka et al., 2015). Comparative 
proteome analysis of drought-responsive proteins during the seedling stage in two cotton (G. 
hirsutum) cultivars, drought-tolerant KK1543 and drought-sensitive Xinluzao26 have performed 
proteomic changes to response for drought in roots of cotton plants. These divergently expressed 
proteins may be related to biochemical pathways that are responsible for drought tolerance 
(Zhang et al., 2016). 

However, it is still required to study in-depth to find the most precise indices to evaluate plant 
drought tolerance. It is important to elucidate and induce drought-tolerant traits via metabolome 
analysis and identify their traits or parameters as potential indices for crop improvement. 

This paper investigated the relationship between metabolites and morphophysiological 
parameters with yield components of cotton cultivars 'Ishonch' and "Tashkent-6' under water 
deficit condition. Results indicate that 'Ishonch' has shown more drought tolerant genotype while 
'Tashkent-6' is drought sensitive. Conclusion of this study is the characterization of cotton 
cultivars by using morphological, physiological and metabolomic approaches can provide 
valuable information for selection criteria in plant breeding programs. 


2 Materials and methods 


2.1 Experimental site, irrigation condition, and plant material 


The experiment was conducted in Zangi-Ata, a district in the Tashkent region of Uzbekistan. The 
Tashkent region has a continental climate, according to the Koppen climate classification system 
(Peel et al., 2007). Hence, the region experiences cold winters and long, hot, and dry summers. 
The Zangi-Ata experimental station is located at 41.178°N and 69.129°E, with an elevation of 
398 m a.s.l. The annual photoperiod (light/dark) is 16/8 h. Temperatures increase in April during 
cotton sowing season and decrease in late September before the harvesting period. Sunny days are 
between 175-185 d. Rainfall varies from 6 to 70 mm in the dry season for a period of 5—6 months 
(Fig. 1). The crops require intensive irrigation throughout this vegetation period. 
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Fig. 1 Monthly average temperature and precipitation in Tashkent region, Uzbekistan. The data include 1981- 
2010 (a) and 2014 (b). Source: Centre for Hydro-meteorological Services, Tashkent. Max. Temp, maximum 
temperature; Min. Temp, minimum temperature. 


An optimal cotton irrigation protocol is widely-used for agriculture in Uzbekistan. Cotton is 
irrigated according to the 1-2-1 (pre-flowering—flowering—boll opening) sequence with 900 
m*/hm? of water applied before flowering, two applications of 1200 m*/hm? each during 
flowering, and 900 m*/hm? prior to boll opening phases. Soil moisture also contributes water 
during seed germination. A modified irrigation protocol was also developed for deficient 
irrigation conditions. It has a 1-1-0 sequence which limits water availability during flowering and 
boll development stages and reduces the total irrigation requirement to 2100 m?/hm? of water. 
This study evaluates the response of selected cotton (G. hirsutum) cultivars, 'Ishonch' and 
"Tashkent-6', to water deficiency under field conditions. Both genotypes have an average fiber 
production (3.5-3.8 tons/hm?) but different tolerances to drought. Cotton cultivars were planted 
20 cm apart in 90 cm wide by 50 m long furrows on 2 April, 2014. Full and deficit irrigation 
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conditions were separated by a specified distance. Insecticides Bi-58 (BASF, Germany) and 
Hexachloran were applied to aphids and cotton bollworm, respectively. Seasonal application of 
fertilizer per annum consisted of 250 kg/hm? nitrogen, 180 kg/hm? phosphorus and 115 kg/hm? 
potassium. Fertilizers were applied during tillage or before irrigation. Morphological, 
physiological, and economic traits of selected cultivars were studied under water deficit 
conditions. Parameters such as plant height, leaf area (LA), weight of dried leaf (WDL), 
chlorophyll fluorescence (CF), relative water content (RWC), excised-leaf water loss (EWL), 
specific leaf weight (SLW), transpiration rate (TR) and metabolite profile were used to monitor 
stress conditions as compared to full irrigated conditions. 


2.2 Measurement of morphological traits 


Plants with the same average height were randomly selected from each individual plot in fifteen 
different locations for morphological measurements. Plant height was determined by using a 
standard meter stick to measure from the ground to the top of the plant. Numbers of internodes, 
sympodial and monopodial branches were manually counted. The mass of boll containing seeds 
and lint was weighed before and after manually separating seeds. Then, lint percentage (lint %) 
and lint index were calculated. Cultivar bolls were considered mature in the plant if more than 50% 
of bolls were open at the time of harvesting. Bolls from the first sympodial branch in the first 
node were selected to analyze boll mass, seed and fiber traits. Fifteen biological replicates of each 
genotype were used for all measurements. 


2.3 Measurement of physiological parameters 


Relative water content (RWC) was measured following the method proposed by Smart and 
Bingham (1974). Fully expanded stem leaves were randomly plucked and immediately placed 
into a container at 10°C—15°C. After 30 min, the RWC was measured. Leaves were weighed 
before and after placing fresh leaves. Leaves were dried for 24 h at 70°C and reweighed to 
calculate the new RWC. Water loss was determined using the Amid et al. (2013) method. 
Transpiration rate was determined by Maniou et al. (2014). Leaf surface density was calculated as 
described in Garnier et al. (2011) and chlorophyll content was measured using a SPAD502 plus 
(Minolta Corp., Osaka, Japan). 


2.4 Non-targeted leaf metabolite profiling 


For metabolite profiles, the 40% methanol extraction procedure was employed as described by 
(Keinänen et al., 2001). Approximately 100 mg of homogenized fine powder leaf tissue was 
placed in a 2-mL Eppendorf tube with 1 mL of extraction buffer (40% MeOH (v/v, 
methanol/water)) and 50 mM acetate buffer (40 mM acetic acid and 44 mM ammonium acetate, 
pH 4.8) and centrifuged at 160,000 g for 10 min at 4°C. The supernatant was transferred into a 
new tube. Totally 15 uL of the extract was transferred into a glass vial with inserts and 2 uL of 
40% methanol leaf extracts were separated using HPLC (Agilent Technologies-1260, USA). An 
Eclipse XDB-C18 narrow-bore analytical column (3.5 um, 2.1 mmx150 mm, Agilent 
Technologies, USA) was used with an applied separation gradient having the following 
parameters: 0 to 0.5 min, 80% A (de-ionized water, 0.1% (v/v) acetonitrile and 0.05% formic 
acid), 18% B (acetonitrile and 0.05% formic acid) with 2% C (isopropanol and 0.05% formic acid); 
0.5 to 15 min to reach 72% B+8% C. The flow rate was 250 mL/h. Eluted compounds were 
detected by a 6420 Triple Quad LC/MS (Agilent Technologies, USA) mass spectrometer equipped 
with an electrospray ionization source in negative ionization mode. Typical instrument settings 
were as follows: capillary voltage (4000 V), dry gas temperature (3000 C) and dry gas flow (3 
L/min). Ions were detected from 15 to 2200 m/z at a repetition rate of 1 Hz. Raw data files from 
the 6420 Triple Quad LC/MS were converted to mzData file format and processed using the 
R-package, XCMS script, and CAMERA (http://bioconductor.org/packages). Peak detection was 
performed using the centWave method with parameter settings: ppm=20, snthresh=10, peak 
width=5 to 20 s. The metaboanalyst online software (http://www.metaboanalyst.ca) was used to 
analyze the XCMS output data file and to perform multivariate statistical analyses after 
normalization using Pareto scaling. The mass spectra of compounds were identified using the 
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commercial mass spectral database, MASSBANK (http:/(www.massbank.jp), and the public 
domain mass spectra library from the Max-Planck-Institute for Plant Physiology, Golm, Germany 
(http://gmd.mpimp-golm.mpg.de). 

2.5 Statistical analysis 


Data analysis was performed using StatView (SAS Institute Inc., Cary, NC, USA) with one-way 
ANOVA followed by a Fisher PLSD post hoc test (P<0.05). 


3 Results 


3.1 Water deficiency altered morphological traits and yield components 


Morphological traits were differently affected among the two cotton cultivars under deficit 
irrigated conditions. Among the morphological traits, the plant height was significantly reduced in 
both cultivars under deficit irrigated condition (Fig. 2al, P<0.05—0.01), but the higher reduction 
was observed in "Tashkent-6'. Likewise, the water stress significantly reduced the numbers of 
internodes in both cultivars which correlated to a reduced number of sympodial branches (Fig. 
2a2, P<0.05—0.001). To the contrary, the number of monopodial branches, either under full or 
deficit irrigated conditions, showed no difference (Fig. 2a4). 
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Fig. 2 Water deficit affected morphological traits of the two cotton varieties (cultivars, 'Ishonch' and 
'Tashkent-6'). Plant height, the number of internodes, the number of sympodial branches, and the number of 
monopodial branches were significantly altered in cotton varieties grown under deficit irrigated conditions (a). 
The water stress influence on the economic traits of the 'Tashkent-6' cultivar (b). The bars indicate standard errors 
for fifteen replicates, while the lowercase letters indicate significant differences determined by one-way ANOVA, 
followed by a Fisher PLSD post hoc test (P<0.05). 
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For the water deficit stress, it can be seen that the common economic traits of both cotton 
cultivars were differently affected as well (Fig. 2b). For 'Tashkent-6' under water deficit 
conditions, the cotton yield (P<0.001; Fig. 2b10) was significantly affected through a reduction 
that was correlated with a lessened number of bolls, up to 42.8% (P<0.001; Fig. 2b1). 
Accordingly, the boll mass (lint and seed) was also significantly reduced in 'Tashkent-6', but not 
at all affected in 'Ishonch' under deficit irrigated conditions (P<0.05; Fig. 2b2). Fiber and seed 
traits were also highly influenced by the water deficit stress. Additionally, there was no difference 
observed for the seed number per boll or fiber length in both cultivars, but there was significant 
reduction observed in the seed mass per boll and one hundred seed mass in 'Tashkent-6' 
(P<0.05—0.01; Figs. 2b6—2b8). 


3.2 Effects of water deficit stress on morphophysiological parameters 


Certain physiological parameters were determined to understand crop tolerance to water stress. It 
was found that while both cultivars maintained a leaf water balance, the parameters differed with 
content. For one, it was observed that LA significantly decreased in 'Ishonch' (P<0.001; Fig. 3a), 
but showed no change for 'Tashkent-6' under deficit irrigated conditions. On the contrary, EWL 
showed a substantial reduction in "Tashkent-6' (P<0.05; Fig. 3b), whereas SLW significantly 
increased in 'Ishonch' compared to "Tashkent-6' during water deficiency stress (P<0.05; Fig. 3c). 
Likewise, LDW significantly decreased in 'Ishonch' (P<0.01; Fig. 3d), while no changes were 
observed for 'Tashkent-6'. Additionally, the transpiration rate was significantly lowered in both 
cultivars 'Ishonch' and "Tashkent-6' by 42.5% and 64.4%, respectively (P<0.01; Fig. 3e). Lastly, 
the chlorophyll content was increased in 'Ishonch', but not in 'Tashkent-6' (P<0.05; Fig. 3f), while 
RWC was uniform for both cultivars during deficit water stress (P<0.18—0.82). 
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Fig. 3 Effect of water deficiency on physiological parameters in cotton cultivars under deficit and irrigated 
conditions. The bars indicate standard errors for five replicates and different lowercase letters show significant 
differences determined by one-way ANOVA, followed by a Fisher PLSD post hoc test (P<0.05). 


3.3 Water deficit stress altered metabolite profiles 


The whole metabolite profiles in both cultivars under full and deficit irrigated conditions were 
further analyzed to study plant response to water stress. Totally, 341 metabolite ions features were 
identified during the initial raw data processing step. Principal component analysis (PCA) 
demonstrated metabolites profiles under full and deficit irrigated conditions that were separated in 
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the 'Ishonch' cultivar, but not in the "Tashkent-6' cultivar (Fig. 4a). The PC1 component shows 
variance between different irrigated conditions, while the PC2 component describes variance 
within the groups. The PCA plot clearly demonstrates that metabolite profiles differed in 'Ishonch' 
grown under full and deficit irrigated conditions compared to "Tashkent-6'. Moreover, clustering 
analysis of leaf metabolites under full and deficit irrigated conditions were grouped separately in 
Ishonch', whereas "Tashkent-6' exhibited mixed clustering of metabolites for full and deficit 
irrigated conditions (Fig. 4b). A heat map of selected metabolite levels in the cotton cultivars 
under different conditions indicated various accumulation levels of metabolites in 'Ishonch' and in 
"Tashkent-6' (Fig. 4c). Moreover, the Venn diagram in Figure 4d depicts the different patterns of 
the logarithmic changes of metabolites for full versus deficit irrigated conditions in the cotton 
cultivars. In the result, 14 and 6 metabolites were up- and down-regulated in ‘Ishonch', whereas 
only 5 and 4 metabolites up- and down-regulated in "Tashkent-6'. And 4 and 2 metabolites were 
up- and down-regulated in both cultivars. Finally, a metabolite annotation illustrates the 
involvement of the various metabolic pathways (Fig. 4e), while another contrasts quantitatively 
the cultivars (Fig. 4f). 
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Fig. 4 Metabolite profiling depicting different metabolic responses of cotton genotypes to water deficiency. 
Principal component analysis after Pareto scaling reveals a clear differentiation of leaf metabolite profiles for 
Ishonch' under full and deficit irrigated conditions but not for 'Tashkent-6' (a). Cluster analysis reveals 
differentiation in full and deficit irrigated conditions of the two cotton cultivars (b). Heat map analysis shows the 
comparison of the logarithmic (log2) fold changes in metabolite abundance in the water stress compared to full 
irrigated conditions in the cotton cultivars. Log2-transformed values of 'Tashkent-6' were sorted from lowest to 
the highest value. Metabolite abundance ratios with two-fold change >2 for up-regulated probes and < —2 
down-regulated (Student f¢ test, P<0.05) (c). A Venn diagram represents the number of log2-fold changed 
metabolite ions in the two cotton genotypes under deficit conditions (d). The sum of predicted compounds 
belonging to metabolic classes (e) and their logarithmic fold change regulation in the two cotton cultivars during 
water deficit stress compared to full irrigated conditions (f). 


Figure 5 depicts the drought response from selected metabolites. Observation of the results 
shows proline and galactinol significantly increased in 'Ishonch' under deficit irrigated conditions, 
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while the levels of putrescine, tyrosine, benzoate, cyanidin significantly increased in both cultivars, 
but differed in accumulation. Additionally, the levels of inositol and glucosinolate were reduced 
only in 'Ishonch' cultivar under water deficit conditions (Fig. 5). 
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Fig. 5 Mass spectrometry data of measured metabolite levels in the cotton cultivars grown under full and deficit 
irrigation conditions. The graph charts show relative levels of selected metabolites in full and deficit irrigated 
conditions. Different lowercase letters indicate significant differences as determined by one-way ANOVA, 
followed by a Fisher PLSD post hoc test (P<0.05). 


3.4 Relationships among physiology, metabolite and yield under deficit irrigated condition 


Relationships among yield components and physiologic and metabolic parameters were examined 
by Pearson’s correlation (Tables 1 and 2). The findings in Table 1 show that yield and the 
morphological traits differed under water stress conditions. The plant yield was highly dependent 
on the number of bolls that had a strong positive correlation under both irrigated conditions 
(Table 1). However, the plant yield was not significantly correlated with the number of bolls 
under either full or deficit irrigated conditions. Similar results were observed between yield and 
the opened bolls. 


Table 1 Correlation coefficients between common economic traits and yield of the cotton cultivars under full 
and deficit irrigated conditions 


Index Full Deficit Full/deficit Deficit/full 
Seed weight per 100 0.912" 0.689* 0.895* 0.570* 
Seed mass 0.872* 0.697* 0.724* 0.491* 
Boll mass 0.928" 0.802" 0.684" 0.487" 
Lint % 0.694° 0.868" 0.701* 0.867" 
Lint mass per boll 0.889" 0.815" 0.727" 0.177 
Total bolls 0.859" 0.944" 0.451 —0.048 
Opened bolls 0.935* 0.972* 0.284 0.078 
Internodes 0.816" 0.173 0.895" 0.383 
Sympodial branches 0.767" —0.189 0.810" 0.014 


Note: * indicates significance at P<0.05 level. 


The relationship between number of bolls and the number of sympodial branches was 
positively correlated under full irrigated conditions but not under deficit irrigated conditions 
(Table 2). A similar relationship was also observed for number of internodes. Seeds, lint traits, 
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and boll mass were positively correlated with yield under both conditions. Water stress highly 
decreased the number of bolls that had an impact on plant yield, but had less influence on the 
correlation between yield and boll traits. Relationships amongst economic traits were strongly 
positive under full irrigation conditions but little or no correlation was seen for internodes and 
sympodia with other traits under deficit irrigated conditions (Tables 1 and 2). 


Table 2 Correlation coefficients between common economic traits and yield of cotton cultivars under full 
(lower diagonal) and deficit (upper diagonal) irrigated conditions 


Y SW SN SM BM L% LIMB FL FI TB © I S PH 

Yield (Y) 0.69* 0.81 0.69" 0.80" 0.87" 0.82" 0.92" 0.85" 0.94" 0.97" 0.17 -0.19 0.85* 
Seed weight (SW) 0.91" 0.90* 0.87" 0.90° 0.92" 0.92* 0.79" 0.89* 0.77* 0.65 0.73* 0.53 0.88" 
Seed number (SN) 0.91” 0.88" 0.94" 0.96" 0.95" 0.97* 0.89* 0.91" 0.88* 0.80* 0.51* 0.30 0.90" 
Seed mass (SM) 0.87" 0.90* 0.86" 0.93" 0.86" 0.93* 0.81* 0.79" 0.80* 0.69° 0.51 0.36 0.81° 
Boll mass (BM) 0.92" 0.92". 0.87" 0.88" 0.93" 0.98" 0.84" 0.86" 0.89" 0.81" 0.44 0.25 0.86" 
Lint % (L%) 0.69" 0.86" 0.63" 0.81" 0.77" 0.96" 0.93" 0.97" 0.92" 0.85" 0.51 0.26 0.96" 
ees 0.88" 0.78" 0.91* 0.73* 0.80° 0.49 0.86" 0.89° 0.89° 0.80° 0.50 0.28 0.90" 
Fiber length (FL) 0.58 0.73" 0.53 0.78" 0.69" 0.94" 0.36 0.92" 0.94" 0.92" 0.34 0.08 0.88" 
Fiber index (FI)  0.91* 0.84" 0.89* 0.77" 0.89" 0.54 0.85" 0.41 0.89" 0.82" 0.57 0.28 0.93" 
Total bolls (TB) 0.85" 0.72" 0.92" 0.74" 0.79" 0.40 0.86" 0.31 0.88" 0.95° 0.25 -0.06 0.87" 
Open bolls (OB)  0.93* 0.80" 0.89" 0.81" 0.84" 0.48 0.87" 0.35 0.94" 0.93° 0.06 -0.23 0.82* 
Internodes (1) 0.81" 0.82" 0.83 0.89 0.82" 0.73" 0.68" 0.71" 0.72" 0.75" 0.75" 0.85" 0.52 
Sympodiya (S) 0.76" 0.71* 0.91* 0.80° 0.75" 0.43 0.77" 0.42 0.80" 0.91* 0.81" 0.77" 0.26 


Plant height (PH) 0.52 0.28 0.61 0.34 0.37 -0.18 0.64 -0.24 0.60 0.80° 0.72" 0.38 0.71 


Note: Each genotype was analyzed as a data point (n=10). * indicates significant correlation at P<0.05 level. 


This study analyzed the correlation between yield and physiological and metabolites parameters 
(Table 3). The relationship between proline and yield in 'Ishonch' under deficit irrigated condition 
revealed positive and moderate correlation, but no correlation was observed for 'Tashkent-6' 
(P<0.05 and P<0.16, respectively). Correlation between putrescine and yield was negative for 
both 'Ishonch' and 'Tashkent-6' (P<0.08 and P<0.14, respectively) under full irrigation. However, 
no correlation was observed for 'Ishonch' under deficit irrigated condition. Galactinol was 
negatively correlated with yield in 'Ishonch' under full irrigated condition (P<0.07). However, the 
relationship between these traits was strongly and positively correlated under deficit irrigated 
condition (P<0.01). No correlation between galactinol and yield was observed in 'Tashkent-6' 
under water deficiency. The anthocyanin content is expressed as cyanidin-3-glucosides belonging 
to phenylpropanoid compounds that commonly accumulated in cotton. The relationship between 
cyanidin-3-glucoside and yield was positively correlated in both cultivars under water stress 
(P<0.23 and P<0.14, respectively). Correlation analysis amongst morphophysiological traits 
revealed various relationships (Table 4). LA together with SLW is one of the main 
morphophysiologic parameters that controls water loss (Xu and Zhou, 2008). LA is negatively 
correlated with SLW under full irrigated condition (P<0.05) and there was no correlation under 
deficit irrigated condition in 'Ishonch' (Table 4). However, these traits in "Tashkent-6' were strong 
negatively correlated under both conditions (P<0.01). A negative relationship was observed 
between LA and water content in both genotypes under full and deficit irrigated conditions 
(P<0.01—0.05). However, LA weakly and negatively correlated with EWL in 'Ishonch' under both 
conditions (P<0.08). Apparently, this correlation in "Tashkent-6' was positive under deficit 
irrigation (Table 4). This indicates that unchanged LA 'Tashkent-6' that assists in leaf water loss. 
Transpiration rate was strongly correlated with EWL in both cultivars under full irrigation 
(P<0.01). However, it was negative in 'Ishonch' and positive in "Tashkent-6' under deficit irrigated 
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condition (P<0.01 and P<0.05, respectively; Table 4). In 'Ishonch’, the correlation between RWC 
and EWL was high positive under full irrigation, and moderate and positive under deficit 
(P<0.01) irrigation. Apparently, these traits in 'Tashkent-6' was moderate and positive under full 
condition and negative under deficit irrigated condition (P<0.01). Correlation analyses revealed 
that yield was correlated with LA, SLW, DLW, RWC, CC, and EWL but not differed in 
coefficients in Ishonch' under both irrigated conditions. However, transpiration intensity was 
positively correlated under deficit irrigated condition. Apparently, "Tashkent-6' resulted in 
different responses in relationships between yield and physiological traits for full and deficit 
irrigated conditions. From this paragraph, we infer that 'Ishonch' compared to "Tashkent-6' is less 
sensitive to water deficiency. 


Table 3 Correlation between yield and physiologic/metabolic attributes of two cotton cultivars under full and 
deficit irrigated conditions 


Yield or Ishonch Tashkent-6 
metabolite traits Full Deficit Full Deficit 
Proline 0.122 0.380 0.092 —0.098 
Putrescine —0.457 —0.108 —0.843* —0.550 
Galactinol —0.494 0.756" 0.641* 0.049 
Inositol 0.478 0.140 0.424 —0.571 
Tyrosine —0.611 0.078 —0.126 0.557 
Benzoic acid 0.041 0.095 0.932* 0.898* 
Cyanidin 0.336 0.570 —0.116 0.468 
Glucosinolate —0.856 0.031 —0.704* 0.578 
Transpiration —0.915* —0.421 —0.724* —1.000* 
Leaf area 0.106 0.518 —0.789* —0.976* 


Note: Each genotype was analyzed as a data point (n=10). * indicates significant correlation at P<0.05 level. 


Table 4 Relationships among morphophysiologic traits in 'Ishonch' (lower diagonal) and ‘Tashkent-6' (upper 
diagonal) under full and irrigated conditions 


Y LA SLW DLW RWC EWL TR CC 
0.70* —0.99* 0.07 —0.73* —0.97* —0.60* 0.86* 
Yield (Y) 
0.99* —0.78* 0.95* —0.98* 0.56 0.37 0.86* 
0.99% —0.81* 0.76* —1.00* —0.50 0.15 0.97" 
Leaf area (LA) i . 
0.99%? —0.86* 0.90* —1.00* 0.67 0.50 0.78* 
—0.62* —0.72* —0.23 0.80* 0.91* 0.46 —0.94* 
Specific leaf weight (SLW) ; 
—0.57 —0.44 —0.56 0.90* —0.95* —0.87* —0.36 
0.85* 0.77* —0.11 —0.70* 0.19 0.76* 0.56 
Dry leaf weight (DLW) i i l 
0.64* 0.74° 0.27 —0.90° 0.29 0.08 0.97” 
Relative water content el tae 0.04 -1.00 0.54 -0.10  —0.98* 
(RWC) -1.00"  —0.98* 0.61 —0.60 —0.72" -0.55 ~0.74" 
Excised-leaf water loss -045 -0.33 -0.42 -0.86* 0.90" 0.78" -071° 
(EWL) -0.55 0.42 1.00* 0.29 0.60" 0.98" 0.06 
—0.61 —0.50 —0.24 0.94" 1.00 0.98" —0.11 
Transpiration rate (TR) ; 
0.67* 0.56 —0.99* —0.15 —0.70* —0.99* —0.15 
Chlorophyll content (CC) —0.27 —0.40 0.92* 0.28 —0.30 —0.73* —0.60 
—0.79* —0.87* —0.05 —0.97* 0.80* —0.07 —0.08 


Note: ° indicates the full irrigated condition; ” indicates the deficit irrigated condition; * indicates significance at P<0.05 level. 
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4 Discussion 


Drought, or water stress, is considered as a limiting factor in agricultural productivity (Boyer, 
1982). The generation of drought-tolerant crop plants is a challenge for plant breeders due to the 
lack of region-specific drought-related indices. Several studies have shown that drought stress 
significantly influences morphology (Pettigrew, 2004; Jaleel et al., 2009). However, 
morphological changes could vary due to different levels of water supply (Basal et al., 2009). In 
this study, the usual watering scheme was modified as described above (see material and 
methods). Change the scheme of traditional irrigation protocol alters crop productivity. However, 
assigning drought indices is important in plant breeding during water deficiency. 

Plant height could serve as an index for drought tolerance with respect to plant morphology, but 
not explicitly for plant breeding. Several studies have proven that water stress highly influences 
plant stems (Pace et al., 1999; Pettigrew, 2004). Recently, McGarry et al. (2016) showed that 
dynamic changes in Single Flower Truss (SFT) and Self-Pruning protein levels guide meristems 
between monopodial and sympodial branching programs in upland cotton (G. hirsutum). Further 
expressions of the genes were up-regulated in sympodial and monopodial stems rather than in 
leaves. Functional analyses of these genes show that they regulated indeterminate growth of 
sympodial and monopodial branches. In all likelihood, drought could alter the gene expression of 
these proteins that results in the change of plant stalk (Fig. 2). Moreover, Pettigrew (2004) 
reported that branching and stem development typically attenuated in drought stress result in 
reduction of bolls. Most of them are retained at the first fruiting position. 

Plant yield is dramatically affected by limited water supplies due to reduction boll number. 
Grimes et al. (1969) demonstrated a strong positive relationship between yield and boll number 
that is consistent with current results. However, water stress changed correlation between these 
traits that could be considered as an important index in order to evaluate drought tolerance. 
Additionally, boll mass, the numbers of seed and lint traits have also proven to influence plant 
yield (Basal et al., 2009). Onder et al. (2009) reported that drought affects seed mass but not seed 
number which is consistent with our results. Evidently, drought stress reduces boll number due to 
the production of fewer flowers (Grimes, 1969; Gerik et al., 1996). However, ‘Ishonch' 
demonstrated non-significant decreases of bolls trait which suggests that flower development was 
not affected by water stress compared to 'Tashkent-6'. 

Many studies have reported that water deficiency can dramatically affect cotton yield during 
boll development (Plaut et al., 1992; Radin et al., 1992; Trolinder et al., 1993), which is 
considered as the most water-deficit-sensitive period for cotton. However, there is still debate 
over the definition of the most sensitive period for the vegetative growth of cotton. Limited water 
supply during pre- or post-flowering stage demonstrated the influence formation of a number of 
bolls which reflects on yield. Our modified watering scheme was designed to one-time irrigation 
during pre-flowering period that showed a significant reduction in common economic traits for 
"Tashkent-6' compared to 'Ishonch'. The result from this study determined that flowering period 
was sensitive to drought. Also, fiber production, which is based on the number of ovules 
contained per boll and flowering stage, was a critical period for yield determination. Studies have 
shown that pollen development and pollen tube growth are sensitive to abiotic stresses due to the 
environment (Burke, 2002; Snider et al., 2011). Consequently, affected flower development, 
during water deficit stress, influences the boll development. Results are consistent with respect to 
the work from McMichael et al. (1973), who demonstrated that young bolls abscise if water stress 
occurs during fourteen days after anthesis. All of these morphological alterations could be caused 
by the changes in water balance of plants. 

The impact of deficit irrigation was determined for a number of interrelated physiological 
parameters which include cellular water content, stomatal conductance, leaf weight, transpiration 
intensity, and photosynthesis (Meyer et al., 2014). High RWC in leaf tissue has been reported as a 
physiologic selection criterion in breeding for drought tolerance (Tahara et al., 1990). However, 
EWL can also be considered as the selection criterion for the characterization of plants for 
drought tolerance. EWL is genetically regulated which can also be applicable in breeding. The 
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level of EWL in 'Ishonch' demonstrated less water loss compared to "Taskent-6'. Dhanda and Sethi 
(1998) demonstrated that EWL is genetically inherited in winter bred wheat. Similarly, 
hybridization analysis of cotton varieties with 'Ishonch' and their offspring demonstrated less 
water loss and indicated genetic inheritance of this trait in second generations. Additionally, 
regulation of leaf area during water stress also plays an important role in water loss. 

Regulation of leaf area controls transpiration intensity (Xu and Zhou, 2008). Alves and Setter 
(2004) reported that leaf expansion under drought stress is decreased due to the reduction of 
transpiration rate. Decreased leaf area in 'Ishonch' might tend to close stomata, or reduce them, in 
order to avoid water loss. Moreover, Pettigrew (2004) reported the reduced leaf area during water 
deficits, but also noted an increase in specific leaf weight (SLW). Specific leaf weight (SLW), or 
leaf thickness, is the integral structure of the leaf which is also known to have a correlation with 
the photosynthetic rate. Evidently, 'Ishonch' increased leaf thickness; however it also reduced leaf 
area when compared to "Tashkent-6'. Higher leaf thickness and reduced leaf area in 'Ishonch' 
could prevent water loss during drought stress periods which, therefore, exhibits drought tolerance. 
Furthermore, ‘Ishonch’ could induce stomata closure during water stress. Patil at al. (2011) 
concluded that higher SLW might have enhanced the photosynthetic rate, thus resulting in better 
growth. 

It is understood that plants change their metabolome during abiotic stresses (Obata and Fernie, 
2012; Witt et al., 2012; Meyer et al., 2014), which suggests that these changes could be used to 
characterize plants for stress tolerance (Duan et al., 2012; Saxena and Cramer, 2013). Drought 
tolerant plants could also reorganize its metabolome profiles during water deficiency by 
producing or increasing the number of secondary metabolites. Metabolites, such as proline, 
galactinol, benzoic acid, inositol, tyrosine, cyaniding, and glucosinolates (Fig. 5) have shown to 
be significantly altered in plants during drought stress period (Taji et al., 2002; Nishizawa et al., 
2008; Duan et al., 2012; Obata et al., 2015). Obata et al. (2015) showed that candidate metabolite 
markers significantly correlated with yield and field trials in maize during drought stress. 
Genetically, 'Ishonch' was originated from a hybrid between line L-27 (G. hirsutum and G 
barbadense) and "Tashkent-6' where L-27 manifest drought tolerance (Axmedjanov et al., 2016). 
However, 'Tashkent-6' was originated from G. hirsutum ssp. mexicanum and C-4227 (G. hirsutum) 
with subsequent backcrossing with C-4227 (Mirakhmedov, 1974). Analysis of different G 
hirsutum and G. barbadense accessions for drought tolerance revealed that G barbadense had a 
greater potential for volatility and acclimatization to unusual conditions, capable of better 
resisting the adverse effects of salt (Tiwari et al., 2013) and drought (Mammadova et al., 2015). 
This could be one of the reasons why 'Ishonch' had less water loss. Correlation between proline 
and yield showed positive relationship under water stress conditions with respect to drought 
tolerant 'Ishonch'. Early work with cotton demonstrated that proline amounts highly accumulated 
during drought stress (McMichael and Elmore, 1977). Proline plays an important role as a 
regulatory mechanism of water loss by reducing cell water potential (De Fatima Fumis and Pedras, 
2002). Leaf water potential was reduced under water stress conditions. However, the level of free 
proline depends on the genetic background of the plant. Guimarães et al. (2008) demonstrated that 
the amount of free proline was dependent on sugarcane genotypes. Moreover, increased level of 
galactinol strongly correlates to yield under stress condition. In addition, galactinol could have a 
novel function in protecting plant cells from oxidative damage caused by stress (Nishizawa et al., 
2008). Over-expression of gene encoding galactinol synthase in Arabidopsis resulted in the 
accumulation of galactinol that subsequently enhanced tolerance to drought (Taji et al., 2002). 
These findings might validate that 'Ishonch' is genetically drought tolerant which is and able to 
re-program its primary metabolism to secondary metabolism. 

Interestingly, the benzoic acid level had a strong positive correlation with yield in "Tashkent-6' 
at both watering schemes. On the contrary, no correlation observed in 'Ishonch', suggesting that 
this phenomenon is genotype specific and hence, has no relation to watering conditions. The 
benzoic acid level could, however, promote an increased sensitivity to oxidative stress which is 
promoted by water deficit. It was shown that salicylic acid (SA) synthesized from phenylalanine 
via benzoic acid which is significantly associated with high levels of H202 in A. thaliana, which 
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resulted in a significant oxidative damage while the SA-deficient transgenic line showed 
enhanced tolerance to salt and osmotic stress in the absence of SA accumulation (Borsani et al., 
2001; Mohamed-Shater et al., 2015). Moreover, the inhibition of benzoic acid can control the 
endogenous response to salt stress and prevent SA accumulation in rice under oxidative stress 
(Sawada et al., 2006). Conversely, the benzoic acid directly associated with plant yield. The 
application of benzoic acid showed a significantly positive effect on fruit yield of different crops 
(Benavides-Mendoza et al., 2012; Abouziena et al., 2016). Our results showed that 'Tashkent-6' 
significantly outperformed 'Ishonch' by yield components (number of bolls, boll mass, lint weight, 
seeds per boll, seed weight, lint index, and cotton yield) under full irrigation scheme. This might 
be attributed to genotype-specific regulation of high benzoic acid level in "Tashkent-6'. 

Under deficit irrigation scheme, including strong water limited stress conditions, it is clearly 
seen that metabolite profiles in 'Ishonch' demonstrated significant separation from the primary 
metabolic state compared to "Tashkent-6'. Results were consistent with earlier work by the study 
group. This study dissected experimented metabolome profiles in different conditions of water 
supply (1-2-1, 1-1-0, 1-0-0, and 0-0-0) (Bozorov et al., 2016). 

Common economic, morpo-physiologic, and metabolic traits responded differently in the 
selection of two cotton cultivars under deficit stress conditions compared to full irrigation. 
Consequently, those traits correlate each other differently as well. Correlation analyses between 
yield and economic traits, physiologic and also with metabolic parameters, demonstrated that 
Ishonch' is more drought tolerant than 'Tashkent-6'. Yield is a vital index in cotton breeding 
which must be either positively or negatively correlated to other physiologic, metabolic, and 
morphologic traits that define right assessment of lines in cotton improvement. Our results 
revealed that in 'Ishonch’, the correlation between yield and morphologic, metabolic, physiologic 
traits were not highly differed but 'Tashkent-6' shows varied correlation coefficients under both 
conditions. Uniformity, correlation coefficients amongst traits could be considered as the stress 
tolerant proxy, which are important indices in plant breeding. 


5 Conclusions 


In this work, we studied morphologic, physiologic, and metabolic responses to water deficiency. 
Current findings demonstrated that decreasing leaf area and increasing leaf thickness prevent 
water loss which depends on genotype. The chlorophyll content tends to increase slightly when 
water stress applied in 'Ishonch'. Specific leaf weight and transpiration rate were significantly 
changed during the water stress in cotton genotypes. Correlation between yield and most of 
physiologic traits was uniform in ‘Ishonch' compared to 'Tashkent-6'. Yield component, and 
number of bolls which highly depend on flower abscission was not significantly changed in 
‘Ishonch'. Hence, in conclusion, the flowering stage is sensitive to drought stress which reduces 
the rate of flower abscission. The theory is promising to consider flowering time, physiologic, and 
metabolic indices in the characterization plant for drought tolerance. Mass spectrometry based 
fingerprinting is reliable for rapid screening cotton genotypes that can be served as a selection 
index for drought tolerance. Moreover, Ishonch' can be used as a donor for drought tolerance 
while "Tashkent-6' as a donor for high productivity in cotton breeding programs. 
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